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Using unimproved and non-perturbatively 0(a) improved Wilson fermions, results are given for the 
three lowest moments of unpolarised nucleon structure functions. Renormalisation, chiral extrapolation 
and the continuum limit of the matrix elements are briefly discussed. The simulations are performed for 
both quenched and two flavours of unquenched fermions. No obvious sign of deviation from linearity in 
the chiral extrapolations are found. (This is most clearly seen in our quenched unimproved data, which 
extends to lighter quark mass.) Possible quenching effects also seem to be small. The lowest moment thus 
remains too large, so it seems to be necessary to reach smaller quark masses in numerical simulations. 
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1. INTRODUCTION 

Much of our knowledge about QCD and the 
structure of the nucleon has been derived from 
Deep Inelastic Scattering (DIS) experiments, 
either IN ^ IX {I = e~,/i~) via the ex- 
change of a photon or vin — > IX ^ Vip — > IX 
via or W~ respectively. The cross sec- 
tion is determined by the structure functions 
F{x, Q^) = 2xFi,F2 and additionally for neu- 
trino DIS, xF^. Considering non-singlet {NS 
or FP—F"^) combinations only to avoid any ex- 
tra gluon terms, the operator product expan- 
sion relates moments of these structure func- 
tions to nucleon matrix elements Vn-NS as 
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where Ep°J .^^ is the Wilson coefficient and 
v^?]^g is proportional to the matrix element of 
_ Of''). 0(«) is a quark bihnear form, in- 
volving a 7-matrix and n — 1 covariant deriva- 
tives (see 1^ for our conventions). This 'renor- 
malisation group invariant' (RGI) form gives 

*Talk given by R. Horsley at LatOl, Berlin, Germany. 



a clean separation between a non-perturbative 
number v^!j^g (which can be computed on the 
lattice) and a. function EfP^^.j^g{Q) (which is 
perturbatively known). This is a possible di- 
rect comparison between the experimental re- 
sult and the lattice result. More practical at 
present, however, is to use parton density func- 
tions (e.g. MRS, 1^) determined from global 
fits and related to the structure function by a 
convolution. Finally, in a scheme S at scale 
AZliM)vlj,giM) with 
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which in the MS scheme with A*^^ from |^ 
at M = 2GeV gives for quenched 0.732(9), 
0.596(10), 0.534(13) (n = 2, 3, 4 respectively) 



while for unquenched, [A.Z, 
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0.695(10). 



2. LATTICE DETAILS 



The euclideanisation and form of the O^'^ 
has been described in ||l|: V2b can be deter- 
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mined with nucleon momentum p = 0, while 
V2a-, V3 and V4 need a three-momentum with 
one non-zero component. By considering the 
NS term the difficult to compute one-quark- 
line-disconnected terms cancel. 

A given operator can mix with three differ- 
ent classes of operators (with the same quan- 
tum numbers under the hypercubic group): 
[A] of higher dimension; [B] the same dimen- 
sion; [C] lower dimension. At present the 
0(a) improvement operators (class [A]) are 
only known for V2- The associated improve- 
ment coefficients are not completely known, 
Q. For W3 and U4 there are additionally two 
mixing operators belonging to class [B] and for 
V4 a further operator in class [C]. 

For V2 we have found that the numerical 
values of the improvement terms are much 
smaller than the operator itself; thus dropping 
them will cause only an insignificant error. 

Renormalisation can be considered in the 
MOM scheme ~ both perturbatively and non- 
perturbatively, ||], by defining Z^'^" from 

[{q(p)\0-'o-'\q(p))] ^ (q(p)|0«°«"|g(p)), 
at = J? . Perturbation theory gives 
Z^^°"(p,a) = l+ 

where we have now computed B^"'^' (csw) for 
a general value of Cgw , • The perturbative 
renormalisation constant for the first mixing 
term for W3 has also been computed; it turns 
out to be very small. Numerically this term is 
also smaller than the matrix clement so we 
shall also drop it. Problems arise with the re- 
maining mixing terms though: the Born term 
(between quark states) vanishes, making their 
determination using only quark states impos- 
sible. At the present stage of development we 
shall simply ignore these extra terms. 

Practically for the operator renormalisation 
we shall use a variant of 'tadpole improved, 
renormalisation group improved, boosted per- 
turbation theory', {TI-RGI-BPT), 
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(where is the number of derivatives in the 
operator, Uq = {^Tj:U°), = gl/u^ andpi is 
the first perturbative coefficient in Uq). do-i in 
the '□' scheme may be found from the known 
do"i" (from O'^'^'). Results for unimproved 
quenched fcrmions are shown in Fig. and 
compared with the non-perturbative method, 
0. Reasonable agreement is seen. 
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Figure 1. Z^^' using TI-RGI-BPT for unimproved 
quenched QCD (dashed lines) at /3 = 6.0. Also shown 
are the results from jnl, filled symbols. (Expected is a 
window where the numbers are independent of .) 

3. RESULTS 

There have recently been suggestions [p|-pT[| 
for the behaviour of Vn close to the chiral limit. 
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1 - rf„(romp,)2 In ^^.J^^Pli.^^^). 



= AZ^''"'Z^°''0 = AZg(a)0(a), 



rather than a linear form: C„(romps)^ + 
We check this using quenched unimproved 
results (because this data extends down to 
lighter quark mass) in Fig. ||. Although not 
conclusive, it would seem that presently lin- 
ear fits are adequate and any possible non- 
linearities can only show up at rather small 
quark mass outside the present range of data. 
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Figure 2 

fermions at /B = 6.0. A linear fit gives 0.362(7). As 
further examples of fits vie take the unquenched d„ — » 
{3g\ + l)/(47rro/,r)^ ~ 0.663. The dashed-dotted line 
is fi^ = 550MeV (fixed), the long-dashed line fi^ ~ 
250MeV (fitted). The filled square is from MRS j^. 

We first consider the continuum limit using 
quenched 0(a) improved fermions (at (3 = 6.0, 
6.2, 6.4) after taking the chiral hmit. In Fig. |^ 
we show the results. The data for the higher 
moments is unfortunately rather noisy. 
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Figure 3. Continuum extrapolation for quenched 



QCD fOrV^^f^g, U^.^g, 

values 0.344(10), 0.119(38), 0.129(36) respectively. 



and vf^'g , 



with continuum 



Finally we consider unquenched results. In 
Fig. ^ we use a fit function ansatz 
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for 0{a) improved fermions with 7 data sets. 
There is no obvious difference to the quenched 
result - which is surprising, While for 

quenched QCD we expect that the sea term is 
suppressed giving too large a quark contribu- 
tion, this is not so for unquenched QCD. 

In conclusion the results presented here 
would seem to indicate that we have to be 
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Figure 4. Chiral extrapolation (upper picture) and 
continuum extrapolation (lower picture) for v^^^g 
using the data sets in giving a result of 0.373(20). 

much closer to the chiral limit in order to be 
able to perceive the partonic properties of the 
nucleon. Further details will be given in 
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